  function yint=mat4init(x)
    yint=[cos(4*x)
           -4*sin(4*x)];
